Near-field thermal radiation transfer controlled by plasmons in graphene 
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It is shown that thermally excited plasmon-polariton modes can strongly mediate, enhance and 
tune the near-field radiation transfer between two closely separated graphene sheets. The depen- 
dence of near-field heat exchange on doping and electron relaxation time is analyzed in the near 
infra-red within the framework of fluctuational electrodynamics. The dominant contribution to 
heat transfer can be controlled to arise from either interband or intraband processes. We pre- 
dict maximum transfer at low doping and for plasmons in two graphene sheets in resonance, with 
orders-of- magnitude enhancement (e.g. 10 2 to 10 3 for separations between 0.1/jm to lOnm) over the 
Stefan-Boltzmann law, known as the far field limit. Strong, tunable, near-field transfer offers the 
promise of an externally controllable thermal switch as well as a novel hybrid graphene-graphene 
thermoelectric/thermophotovoltaic energy conversion platform. 



Heat transfer between two bodies can be greatly en- 
hanced in the near field, i.e. by bringing their surfaces 
close together to allow tunnelling of evanescent photon 
modes. For two parallel, semi-infinite, dielectric sur- 
faces of index of refraction n, maximum flux enhance- 
ment is known to be n 2 times the Planck's black body 
limiti. However, particularly interesting near-field radi- 
ation transfer phenomena involve thermal excitation of 
various surface modes. Due to their localization and 
evanescent nature, it is only at sub-micron separations 
that these modes become relevant. Measuring near- 
field transfer has been experimentally difficult 2 ^—; nev- 
ertheless, the promise of order-of-magnitude enhance- 
ment over the far field Planck's black body limit has 
made near-field transfer the topic of much research^. A 
promising class of materials for enhancing the near-field 
transfer are plasmonic materials, due to high density of 
modes around the frequency of plasmons. The potential 
of graphene^ as a versatile and tunable plasmonic mate- 
rial has already been recognized in applications such as 
teraherz optoelectronics and transformation optics^—. 
Unlike in metals, where high plasma frequencies make 
thermal excitation of surface modes difficult, plasmon 
frequencies in graphene can be anywhere from the ter- 
aherz to the near infra-red^. In addition, the depen- 
dence of graphene conductivity on chemical potential, 
which in turn can be controlled by doping or by gat- 
ing, allows for a tunable plasmonic dispersion relation. 
Transfer between graphene and amorphous Si0 2 14 i 15 as 
well as application of graphene as a thermal emitter in 
a near-field thermophotovoltaic (TPV) system has been 
reported^. Here we analyze the contribution of plasmon- 
polaritons to graphene-graphene near-field heat transfer. 
The choice of identical coupled systems is predicated on 
the idea that resonant enhancement could lead to even 
greater heat transfer capacity. Indeed, we find maximal 
transfer for resonantly coupled plasmon modes (corre- 
sponding to similar doping in the two graphene sheets), 
which can be orders of magnitude larger than the heat 



transfer between two black bodies in the far field. 

In general, the radiative heat transfer between two 
bodies at temperatures T\ and T 2 is given by 
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H= / cfcj[e(cj,T 1 )-e(cj,T 2 )]/(a;;T 1 ,T 2 ) (1) 
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where 8(w,T) = Tiuj / (e Tlu> / kbT — 1) is the average en- 
ergy of a photon at frequency uj (the Boltzmann factor), 
and /(w;Ti,T2) is the spectral transfer function, charac- 
terizing frequency dependence of the heat exchange (i.e. 
how much heat is exchanged at a given frequency). In 
the context of fluctuational electrodynamics 1 ^, the spec- 
tral transfer function /(w;T!,T 2 ) is calculated in the 
following way: thermal fluctuations in the first (emit- 
ter) medium induce correlations between electric cur- 
rents, which are proportional to the real part of the 
medium conductivity^; next, using Green functions, we 
can find the electromagnetic fields in the second (ab- 
sorber) medium induced by the fluctuating currents in 
the first 1 ^; finally, the radiation transfer is obtained by 
calculating the Poynting flux around (or the ohmic losses 
within) the second medium. This approach has been used 
to numerically calculate the near-field transfer between 
two half-space s 17 ' 20 , as well as generalizations such as two 
slabs^i, sphere and a plane^ 2 - two spheres 2 ^, as well as 
ID periodic structures 2 ^. 

The system we analyze, shown in Fig. [TJ consists of a 
suspended graphene sheet at temperature T\ emitting to 
another suspended graphene sheet held at room tempera- 
ture T 2 = 300 A', and a distance D away. In general, the 
p-polarization spectral transfer function for evanescent 
modes between two bodies is 
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where 7 = \Jui 2 j c 2 — q 2 is the perpendicular wave- 
vector and ri(2) is the reflection coefficient for the bot- 
tom(top) body; note that ri. 2 depend on T, and hence 
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FIG. 1. (a) Schematic diagram of the radiation transfer 
problem: a suspended sheet of graphene at temperature T\ 
is radiating to another suspended graphene sheet at temper- 
ature T2 and distance D away, k-vector components are (7,7, 
for the parallel and perpendicular component, respectively, 
(b) Real and imaginary parts of graphene p-polarization re- 
flection coefficient for fj, = 0.5eV, T = 300K, r = 10~ 13 s. 
Dashed line is the vacuum plasmon dispersion relation Q for 
the graphene sheet. Insets show the real and imaginary part 
of reflectivity at q fts 30eV/hc as a function of uj. 



the T-dcpcndence of f(ui, T±, T2). Integration is over 
the parallel wave- vector q, limited only to the evanes- 
cent (q > uj/c) modes. The spectral transfer function 
@ was derived for the case of two semi-infinite slabs^; 
however, it can be shown that the same expression is 
valid when any of the two bodies is a 2D system, such 
as graphene^. Since graphene absorbs poorly (2.3%) in 
the far field (hence is also a poor emitter) , not including 
the propagating modes is a good approximation. The 
contribution of evanescent s-polarized modes can also be 
calculated using Eq. ^ , but it turns out to be negligible 
compared to p-polarized modes, as we discuss later. We 
assume graphene is completely characterized by its com- 
plex optical conductivity a = ay + icti , which depends on 
angular frequency uj, electron scattering lifetime r, chem- 
ical potential /i, and temperature T . Furthermore, the 
graphene conductivity is taken to be independent of the 
parallel wave- vector q (see discussion below), and consists 
of the Drude (intraband) and interband conductivity, ex- 
pressed respectively as2£ 
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where G(£) = smh(£/k b T)/{cosh((j,/k b T) + cosh{£/k b T)), 
and fj, is the chemical potential. Various electron scat- 
tering processes are taken into account through the re- 



laxation time r. From DC mobility measurements in 
graphene, one obtains^ an order-of-magnitude value of 
t w 10" 13 s. 

First we discuss the electrodynamic properties of a 
single suspended sheet of graphene, inherent in the p- 
polarization reflection coefficient, which is illustrated in 
Figure [T|d. The reflection coefficient is r p = (1 — e)/e, 
where e = 1 + ja/(2eoLo) is the dielectric function of 
graphene^. Its pole e = corresponds to the dispersion 
relation of p-polarized plasmon modest 
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which is shown as the dashed line in Fig. QJ>. Figure Q] 
shows plasmons exist in a strongly non-retarded regime 
(q ^> w/c) , indicating a tightly confined plasmon polari- 
ton mode. Graphene also supports s-polarized surface 
modes with a dispersion relation very close to the light 
line2&. However, due to the large density of states and the 
tightly confined nature of p-polarized surface modes, it 
is the p-polarization that dominates (as our calculations 
confirm) the near-field transfer. 

When two parallel graphene sheets are sufficiently 
close (see Fig. QJi), their plasmonic modes can become 
coupled. The dispersion of these coupled modes is 
1— r\(uj)r2 {oj)e~ 2qD = 0, when q 3> w/c, so 7 w iq, which 
is exactly the pole of the integrand of the spectral trans- 
fer function ([2]). The integrand is illustrated in Figure [2] 
for different values of chemical potential. The coupling 
of modes is strongest when both graphene sheets have 
identical parameters (middle panel in Fig. [2]). In that 
case, their individual dispersions are identical. Never- 
theless, the dispersion of the combined system shows two 
branches that dominate the near-field spectral transfer, 
i.e. the implicit equation 1 — r(uj) 2 e~ 2qD = for to(q) 
has two explicit solutions: Weven(<?) and w dd(9) for the 
even, and the odd mode, respectively. The splitting of 
two superimposed resonances is particularly noticeable at 
smaller wave vectors q. For larger q, the splitting disap- 
pears, and the resonant matching of peaks of Im^ri^) sig- 
nificantly enhances the near-field transfer. As the chemi- 
cal potential of one of the sheets changes (top and bottom 
panel in Fig. [2]) , the plasmons in the two sheets move out 
of resonance, coupling decreases, the peaks in the inte- 
grand approach the individual (vacuum) plasmons dis- 
persion curves, and the heat transfer becomes lower than 
in resonance. 

Figure [3^ shows a highly tunable spectral transfer 
function f p for different values of chemical potential and 
relaxation time. Given the chemical potential, the re- 
laxation time determines which processes (interband or 
intraband) are responsible for the peaks in spectral trans- 
fer. Since interband processes are dominant at high fre- 
quencies, all t curves converge in the high frequency 
limit, where Drude losses are negligible. However, in- 
terband processes can play a leading role even below the 
absorption threshold w w 2/i, particularly for small chem- 
ical potential where thermal broadening of the interband 



3 




0.3 0.5 

u) [eV/h] 



FIG. 2. Contour plot of the integrand (a.u.) in from 
([2]), for two graphene sheets at Ti,2 = 300K, separated by 
D — lOnm. Chemical potentials are fii = 0.5eV, while ^12 is 
different for each plot. Dashed lines correspond to the vacuum 
plasmon dispersion relations for the bottom (1) and the top 
(2) graphene sheet. 
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FIG. 3. (a) Spectral transfer function f p (oj) from for 
plasmons in two graphene sheets at resonance, /tii,2 = A*i 
n,2 = t; Ti,2 = 300if, D = lOnm. Solid green line cor- 
responds to the /i 1(2) = 0.3(0.5)6^, r 1(2) = 10" 13 (10~ 14 )s 
case, (b) Contour plot of the integrated ratio of the near- 
field transfer between two graphene sheets, HgJ, and the far 
field transfer between two black bodies, Hj/ B for plasmons in 
resonance (left, /xi,2 = O.leV) and out of resonance (right, 
/i 1(2) = 0.1(0.3)eV). Here, T 2 = 30QK,ti >2 = l(T 13 s. 



threshold (on the order of few k^T) becomes more signifi- 
cant. For example, for ^,1,2 = O.leV (first peak in Fig.[3ji) 
the similarity between r = 10~ 12 s and r = 10~ 13 s spec- 
tral transfer functions indicates that the majority of loss 
in graphene comes from intcrband processes. On the 
other hand, the Drude (intraband) loss term, usually 
important for oj < fj,, can become dominant at higher 
frequencies, for large enough /1 (third peak). Finally, a 
combination of two loss processes, ^1(2) = 0.3(0. 5)eV, 
T i(2) = 10~ 13 (10~ 14 )s can lead to a hybrid spectral 
transfer. While the use of ^-independent expression for 
graphene conductivity Eq. ^ for intraband processes is 
a good approximation^, one must take care when ap- 
plying Eq. ([3]) to intcrband transitions. As indicated in 
Fig. [JJi, intcrband transitions can play a significant role 
in near-field transfer at low doping levels. Here, the con- 
tribution from the non-zero wave-vector becomes impor- 
tant since it broadens the intcrband threshold from 2/z 
to ~ 2/i — fiqvp. On the other hand, this is similar to 
non-zero temperature effects which also broaden the in- 
tcrband threshold, so we do not expect a qualitatively 
different result with g-dependent conductivity. 

We quantify the heat exchange in the near-field by 
plotting (Fig. [3Jd) the integrated transfer H from Eq. ([T]) 
normalized to the transfer between two black bodies in 
the far field. Factoring in the temperature dependence 
shifts the majority of the near-field transfer to lower 



frequencies, due to the exponentially decaying Boltz- 
mann factor. This implies that while doping or gat- 
ing might be advantageous in some applications (for ex- 
ample, cmitter-PV cell bandgap frequency matching in 
near- field TPV systems^), near-field transfer between 
two graphene sheets is maximized for small values of 
doping, despite the stronger peak in spectral transfer 
for nx,2 = O.ieV vs. fjL 1>2 = O.leV (Fig. 0k). For 
plasmons in resonance with ju 1: 2 = O.leV (left panel, 
Fig. [3)3), we observe orders-of-magnitude increase in heat 
exchange, particularly at small separations (xlOOO for 
D = 20nm, T\ = 800-ftT), but also at separations as large 
as 0.1//ra. At larger separations, we observe (not shown) 
the shift of the peak of the spectral transfer function 
f p to /xi^ = O.leV case (red line in Fig. [3£i), indicat- 
ing that the coupling between highly localized, large q, 
modes becomes weaker, and the transfer is dominated by 
lower-frequency, less evanescent modes. The heat trans- 
fer depends in a complex fashion on the parameters of 
the system, and does not seem to yield a simple func- 
tional dependence on the emitter and absorber temper- 
atures (as is the case for two black bodies). Neverthe- 
less, there is a relative advantage (Fig. ^p) to operating 
at lower temperatures, as the temperature dependence 
of the near-field transfer appears to grow slower than 
the T 4 black body dependence. Finally, we note that 
the temperature dependence of conductivity reduces the 
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resonant effect when two graphene sheets arc at differ- 
ent temperatures. This reduction is more pronounced 
for large temperature difference, shifting the peak of the 
spectral transfer on the order of fc&T; however, the rela- 
tive reduction of the integrated spectral transfer function 
is small, with the main temperature dependence coming 
from the Boltzmann factor. This efficient heat exchange 
between two graphene sheets in the near-field, together 
with recently reported advances in hot carrier extraction 
from graphene^, may offer a potential for a novel, hy- 
brid thermophotovoltaic/thermoelectric solid-state heat- 
to-electricity conversion platform. In addition, this ma- 
terial system could pave the way toward an externally 
controllable thermal switch behavior, where one can, by 



means of doping or gating, tune the resonant coupling 
between the hot and the cold side. 
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